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We calculate the number density of helicity ±3/2 gravitinos produced out of the vacuum by the 
non-static gravitational field in a generic inflation scenario. We compare it to the number density 
of gravitinos produced in particle interactions during reheating. 
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; I. INTRODUCTION 

a^ ■ 

It was realized early on ^ that the gravitino could pose a serious cosniological problem in the context of a hot 
' Big-Bang, if it were once in thermal equilibrium. An unstable gravitino, for instance, would decay in the post-Big- 
bJO Bang nucleosynthesis era, if its mass 7713/2 ^ 10'* GeV, and the entropy produced would ruin the successes of Big-Bang 
^ nucleosynthesis. Similarly, if the gravitino were stable, as e.g., in gauge mediated supersymmetry breaking, its energy 
density would eventually overdose the Universe, if its mass TO3/2 ^ 2keV. A solution to this problem was brought 
. forward in Ref. if inflation took place, the gravitinos present at the time of Big-Bang nucleosynthesis were created 
I ' during reheating, in an abundance possibly much smaller than that corresponding to thermal equilibrium. Cosmo- 
logical constraints on their abundance could then be turned into useful upper limits on the reheating temperature 
^ ' Tr §, typically Tr < 10^ - lO*'' GeV, for an unstable gravitino with TO3/2 < 3 x 10^ GeV, or Tr < 10^ - 10*'' GeV, 
. for a stable gravitino with 1 keV < "I3/2 J5 1 GeV. 

These studies assume that the gravitino abundance has been exponentially suppressed during inflation, and that 
gravitinos were only created in particle interactions during reheating. However, particles can be produced out of 
the vacuum in a non-static gravitational background if their coupling to the gravitational field is not conformal . 
\ A well-known case is the production of gravitational waves or scalar density perturbations during inflation. As we 
t argue below, a massive gravitino is not conformally invariant in a Friedmann-Robertson- Walker (FRW) background, 
0^ ' and our main objective is thus to quantify the number density of gravitinos that can be produced gravitationally 
, during inflation. We will consider a generic inflation scenario within N = \ supergravity and briefly discuss the 
O 1' more particular case of pre-Big Bang cosmology j^]. The cosmological consequences of gravitino production during 
I inflation were briefly discussed in Ref. |Q. These authors did not actually study the gravitational production of 
gravitinos, and rather focused on the spin and spin 1/2 cases, as they were interested in the problem of moduli 
(-H [ and modulini fields. With regards to the gravitino, they assumed that one particle would be produced per quantum 
• • ■ state for modes with comoving wavenumber k Hi, where Hi denotes the Hubble scale of infiation. This estimate 
. ^ ] showed that gravitational production of gravitinos could pose a cosmological problem if the energy scale at which 
■ infiation takes place saturates its observational upper bound, i.e. V^^^ ^ 10*^ GeV, and in this respect, it justifies 
5-H ' further the present work. Our study is more specialized than that of Ref. |^] , as we concentrate exclusively on the 
gravitino. However, it is also more systematic and more detailed, as we derive and solve the gravitino field equation 
to provide quantitative estimates of the number density of gravitinos produced. We also examine different cases for 
the magnitude and dynamics of the gravitino effective mass term during and after infiation. Finally, we also study 
the effect of a finite duration of the transition between inflation and reheating, using numerical integration of the field 
equation. This effect is important, as this timescale defines the "degree of adiabaticity" of the transition, and indeed 
the number density of gravitinos produced is found to be inversely proportional to it. 

The study of the conformal behavior of the gravitino also bears interest of its own, apart from any application to 
cosmology, and to our knowledge, the quantization of quantum fields in curved space-time has been examined for spins 
0, 1/2, 1 and 2, but not 3/2 [Q (although the case of a massless gravitino in a perfect fiuid cosmology was studied 
in Ref. In the present work, we focus on the helicity 3/2 modes of the gravitino. The field equations and the 
quantization of the helicity 1/2 modes are indeed more delicate, due to the presence of constraints. These constraints 
vanish identically if supersymmetry is unbroken |^ ; in broken supersymmetry, these constraints do not vanish, but 
do not induce any inconsistency jlOj. As shown below, these constraints apply to the modes of helicity 1/2, not to 
those of helicity 3/2, and for this reason, we leave the problem of the helicity 1/2 modes open for a further study; 
nonetheless, we present these field equations and their constraints. The number density of helicity 3/2 gravitinos 
produced during inflation that we derive in this paper should thus be interpreted as a lower limit. Quite probably, 
however, the number density of helicity 1/2 gravitinos should be of the same order as that of helicity 3/2, and the 
results correct within a factor of order 2. 
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This paper is organized as follows. In Section II, we derive the gravitino field equation for the ±3/2 helicity 
modes, and in Section III, we calculate the number density of gravitinos produced in a generic inflation scenario. We 
summarize our conclusions and briefly discuss the case of pre-Big-Bang string cosmology in Section IV. All throughout 
this paper, we use natural units H = c = mpi = 1, where mpi = (SvrG)"^/^ is the reduced Planck mass. We note 
Mpi = \/87rTOpi the Planck mass. Furthermore, we restrict ourselves to a FRW background, whose metric is written as: 
ds^ ~ gfj^jydx'^dx'^ = a'^{rj){~dri'^ + dx^ + dy^ + dz^), where a{ri) is the scale factor, and rj denotes conformal time; the 
Minkowski metric is written r]ab- We also use standard conventions on the derivative of the Kahler potential G{zi, z*) 
with respect to the scalar components Zi of chiral superfields: = dG/dzi, Gi* = dG/dz" . Other notations, relative 
to the Dirac matrices, are given in the Appendix. 



II. FIELD EQUATION 

We consider the gravitino in a background of a classical FRW spacetime, in the context of = 1 supergravity, and 
adopt the following lagrangian density: 



In this equation, e represents the determinant of the vierbein e^, R denotes the Ricci scalar, the gravitino 
field, and £,„ represents external matter, more specifically the scalar fields whose dynamics drive the evolution of the 
background metric: we neglect the matter gauge and fermion fields. The gravitino covariant derivative Dp is defined 
as: 

Dp = dp + - (2) 

where ti;°''(e) is the torsion tensor, in which wc do not include ^ torsion, since we neglect the backreaction of 
the gravitino on the metric. Wc included in this covariant derivative the Kahler connection Ap = K^dpZi — Ki-^dpZ*, 
where K denotes the Kahler function. The gravitino is also coupled to matter through the Kahler potential G{z, z*) — 
K{z, z*) +ln (|VF(z)p), with W the superpotential. This term gives rise to an effective mass for the gravitino, which 
we write as m: m = e^^'^. The gravitino field equation can be written in the compact notation [p^ : 

= eM-P-^5^^p^$^ ^ 0, (3) 



where Vp = Dp + ^m'jp. As is well-known [pO|p|, a consistency condition can be obtained by taking the divergence of 
Eq. @, VpR^—Q, which leads, after some manipulations, to: 

[im'-l'' ~ G^'p^^^ + 29^771^^" + 2m75A^a'^'^] = 0, (4) 

where Gp,v is the Einstein tensor, symmetric in the absence of torsion. In this equation, we did not include a term of 
the form e^^ ' d^jXn . . ., since it vanishes in a homogeneous and isotropic background. 

We now define: TZp = Rp — ^jpj'^R^, and rewrite the field equation Eq.(^), in an equivalent way, as TZp — 0: 

7^o = (-rd, +m+ ^7" + ^-f5l°Xo] *o - (do - +n + 7"*. - 0, (5a) 



(5b) 



where H = a' /a. Note that in a FRW background, Xi — 0, and djXp = 0. Another integrability condition can be 
obtained from the difference 'y'^TZq — j'^TZi : 

g'^di'i/j = (fdi -m + 7^7°)7J*j. (6) 

Equation (||) and the constraints Eqs. (Q) and (^ form the system of field equations for the gravitino. 

We now perform a standard decomposition of the gravitino field operator. We rescale the gravitino field, and write: 
^^(3^) = fl('7)~'^^^e)j^c('7, A;)e*'"^; we recall that we reserve latin indices a, 6, c . . ., or a hat, if confusion could arise, for 
Lorentz indices, and = a{r])S'ij^. Note also that ^p{x) transforms with a conformal weight 1/2, and ^a{x) transforms 
with a conformal weight 3/2. Then, we decompose the spatial part of *I'c(?7, k) into helicity eigenstates [O]: 
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*c('7,fc)= C'i,i/2(m + ^;TO,^)e™(fe)V'„,(77,fc), (7) 

m— L, + ,— 

where Ci i/2('7i+ §;m, |) is a Clebsch-Gordan coefficient, and e'", and e~ are polarization vectors. They satisfy: 
gS*gS _ Sss', with s = L, +,— ; in particular, is parallel to k, e+ and e are tranverse to fe, and e+* = e~. 
Similarly, the spinors 4'ms are helicity eigenstates of the helicity operator diag (e'^cr, e'^cr) . More specifically, each 
spinor ipjyis{r],k) is written in terms of a Weyl spinor Xs{k) of helicity s/2, i.e. such that e^crxs = sXsj and mode 
functions h,^g{r], k) and gmsiVi where fc = |A;|, following the notations of the Appendix. In this decomposition, the 

vector-spinors e'"-0L+, e'^V'L-, e^V'H j £""0 hi and form the helicity ±1/2 components of while and 

e~ip are the helicity ±3/2 components. 

The field equation for the helicity ±3/2 modes of the gravitino can now be extracted from Eq. (||). To start with, 
one notes that the helicity ±3/2 components do not appear in the product 7^', because project out the modes 
with spinor helicity ±: e^crx± — 0, and e^crx^ = \/2x±- Therefore, Eqs. (jj), ( |5a| ) and only concern the helicity 
1/2 components, not the helicity 3/2. The field equation for the ±3/2 helicity modes is then obtained by contracting 
Eq. (^) with e^je^r]'"^. This contraction projects out all terms of helicity ±1/2, because these are either parallel to 
k, or of the form 7c . . ., and e^j.e^'y = e^e^ = 0. Thus the helicity ±3/2 components do not mix with the ±1/2 
helicity modes, and their field equation reads: 



^7°ao + nfc + am + ^757° Ao^ ?Ass =0, s = ±. 



(8) 



Finally, this equation can be rewritten in the usual way as two systems of two linear and coupled differential 

equations in the mode functions ft.++(77), g++('7), and h {rj), g (77). For m = 0, and zero Kahler connection, it is 

easy to see that Eq. (H) is identical to the field equation for a massless gravitino in Minkowski spacetime, or, in other 
words, a massless and uncoupled helicity 3/2 gravitino is conformally invariant. 

The gravitino field can now be quantized, following the methods developed for spin 1/2 fermions in curved space- 
time [p^P], or, what is similar, for electrons in an external electromagnetic field |l^. Introducing the shorthand 
notation: ^c±3/2 = e^V'iii it can be checked that the inner product '^\^'q°'^'^i,s, where s = ±3/2, is conserved by 
virtue of the field equations. The solutions of Eq. are normalized according to: V'L^s's' — 6ss' , s, s' — ±, and one 
obtains at all times: 

i'i^T^^'^ts' ^S,,,, (9a) 

¥,MkK'^fAv,-k) = 0, s,s'-±3/2, (9b) 
where the superscript C denotes charge conjugation. The helicity— 3/2 gravitino field operator is written as: 

= / E [fc(fc)*a.(^,fc)e''^^+6t(fc)«'C^(ry,fc)e-'=-], (10) 

^ s=±3/2 

where the 5, 5^ are annihilation and creation operators respectively. They are related by hermitian conjugation as the 
gravitino is a Majorana fermion. Finally, one can relate field operators k) and \['°"'(77, fc), that are solutions of 

the field equation, and whose boundary conditions are respectively defined at conformal times ryin and Jjout, by means 
of a Bogoliubov transform Q| : 

i>T:(,r^,k) ^ aksnMv.k) + (3ks^'^^{Tl~k). (11) 

The Bogoliubov coefficients aks and (3ks satisfy at all times: la^sP + |/3fcsP = 1, as required for a half-integer spin 
field. The occupation number operator for the in quantum state with momentum fc and helicity s, s = ±3/2, in the 
out vacuum, is then |/3/csP, and: 

- \Ks{v)gT{v) ~ Cs{r))hT{v)\^ ■ (12) 

In the following, we solve the field equations for the helicity ±3/2, and use Eq. ( p^ ) to calculate the number density 
of gravitinos produced. 
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III. GRAVITATIONAL PRODUCTION OF SPIN 3/2 



We now assume that the background undergoes an era of inflation, followed by radiation or matter domination. 
The magnitude and the evolution of the gravitino mass term in both epochs are model-dependent, since the scalar 
potential V{z, z*) is tied to the Kahler potential in a non-trivial way: 

V{z, z*) = [Gi. (G^i)}. -3]+D- terms (13) 

Nevertheless, it is well-known that scalar fields generically receive a contribution to their mass of order of the 
Hubble constant |p^ , and we adopt this as an ansatz for the gravitino mass term, i.e. 171 = ^111 during inflation, and 
m = /i2-ff during radiation/matter domination, where fii and ^2 a-re constant parameters, H is the Hubble constant. 
During inflation, H = Hi is also assumed constant. Note that this ansatz may be realised rather generically in infla- 
tionary scenarios. For instance |l5[| , the superpotential y/X(f>'^ gives a potential A^"* (albeit in a global supersymmetry 
approximation), a gravitino mass ~ ^/\(/)^/mp^ (also neglecting Kahler terms), and a Hubble constant ~ ^/XcjP' / m-p\. 
In this model of chaotic inflation, (j) ^ Mpi towards the end of slow-roll, and therefore m ~ H. Similarly, for new 
inflation types of models, the superpotential M^(mpi — cj))'^ /mp\ gives a scalar potential ^ Af* when cf) <C mpi, a 
gravitino mass ~ M^/mpi, and a Hubble constant ~ IvP /mp\. 

The quantities /ii and above can take any value, and, presumably, /^i ^ 1 and 1 HJl^- Note that, strictly 
speaking, this ansatz is justifled as long as /xi_2-ff > ^3/21 where denotes the mass of the gravitino in the true vac- 
uum of broken supersymmetry; provided inflation takes place at an energy scale V^^^ > 10^^ GeV (7713/2/10^ GeV)^/^, 
this relation should be satisfied for reasonable values of /ii and /i2. If, however, V^^^ <C 10^^ GeV (7713/2/10^ GeV)^/^, 
then according to the adiabatic theorem ]l7|,^, the production of gravitinos will be exponentially suppressed. Never- 
theless, for the sake of completeness, we also present results for this case where m is constant during both inflation 
and radiation/matter domination. 

For reasons that are similar to the above, one cannot write a generic Kahler connection Ap for a generic inflation 
scenario. It has actually been argued that if inflation is to proceed via the F— terms, the Kahler function K should 
not have a minimal form [p|Jl^. Out of simplicity, we thus assume that this term is zero. This is realized, for instance, 
in scenarios in which the dynamical scalar field is real. Moreover, as we argue in Section IV in the case of string 
cosmology, a non-zero Kahler connection in a homogeneous and isotropic background does not induce particle creation 
by itself. With these assumptions, the differential equations satisfied by the mode functions read: 

g'L - Wss + {k^ + a^m^ - iskU) g^s = (14a) 

IS 

hss^ {g'ss+iskgss) , s = ±, (14b) 

am 

where 7i has been redefined as 7Y = [am)' / {ain) . These equations can be solved in terms of Whittaker functions 
W\^i^_,-aj{zj) and PF-a.^^^q^ (-z^), j = 1,2, where Zj = 2ikaj\rii\{l -I- (1 -I- fy/lmD/^j); ^ = ±1/2, and rji denotes the 
conformal time of exit of inflation: rji = —H^^, as we set a (77 = 771) = 1. The subscript j — 1, 2 correspond to the two 
eras, j = 1 for inflation, and j = 2 for radiation/matter domination; aj is defined by: 0(77) = (1 + (1 + ri/\rii\)/aj)°'^ , 
i.e. ai = — 1 corresponding to de Sitter, and a2 = 1,2 corresponding respectively to radiation or matter domination. 
The in solution is defined as that which reduces to positive energy plane waves as 77 — > —00, and the out solution as 
that which reduces to positive energy plane waves as 77 — > +00. Using the large argument limit of Whittaker functions, 
one obtains |p^ : 

9++{v) = -^Wl/2,^^,^a,izJ), (15a) 
h'tliv) = ^W^-i/2,.;.,a,(^,), J = 1,2. (15b) 



In this equation j = 1 corresponds to the in solution for 77 < 7/1, and j — 2 corresponds to the out solution. In 
the radiation or matter domination region, the in solution reads g++ = ciiy_i/2,i;^2a2 ("-^2) + £2^1/2,^/^202(^2), and 
the coefficients ci and C2 can be obtained by matching g++ and with g++(?7 < Tyi) and /i++(77 < 7?i) in Eq. ( |T5| ) 
continuously at r/ — jji. Finally, the solutions of helicity —3/2 are expressed in terms of the solutions of helicity 



-1-3/2: h = and g = Using Eq. (12), the asymptotic number of particles produced per quantum state 

\f3ks{il +oo)P then reads: 

|/3fcs(+00)P = 1 \^J-2a2W+l/2,i^,^{Zl)W^l/24^,.,a2iz2) + ^^lW^l/24^,l{Zl)W+l/2,i^,2a2{z2)f , (16) 

i^in^2 
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where zi = —2ik\r]i\, and Z2 = 2ika2\r]i\, and 012 = 1,2 depending on whether inflation is followed by radiation or 
matter domination. The limits /zi ^ 0, or /i2 — > are non-singular, and reduce to the solutions one would obtain in 
either of these limits, even though Eq. (|lj) reads differently in these limits (it decouples into two first-order uncoupled 
differential equations). 




FIG. 1. Plot of the power spectrum {k/Hi)'^\l3ksf of the number density of gravitinos per logarithmic wavenumber interval, 
versus k/Hj, for different transition durations Ar; (771 = —Hf^), as indicated, assuming that ^1 = and fj,2 ~ 1, and that 
matter domination follows inflation. The dashed line, which corresponds to Ar; — 0, is obtained from the analytical solution 
in Eq. (p^; the other curves are obtained from a numerical integration of the field equation, with a smooth transition for 0(77) 
and m{ri) between inflation and matter domination. 

The quantity of direct interest to us is I3/2, defined as the ratio of the number density of helicity— 3/2 gravitinos 
to the entropy density at the time of reheating 77^: 



'\(3kA^ (17) 



where g^s is the effective number of degrees of freedom, g^s — 229 for the minimal supersymmetric standard model, 
Tr is the reheating temperature, or = a(77R), ai = a{r]i), and k = kj Hi — k\'qi\ 



Since |/3fcsP < 1 is imposed by Pauli blocking, the integral in Eq. (17) is dominated by the high wavenumber modes. 
As a matter of fact, this integral diverges linearly, since |/3fesP oc fc"'^ for fc|7/i| ^ 1, according to Eq. (p^. This 
divergence is unphysical, and results from the sudden transition approximation, as is well-known, for instance, in the 
case of gravitational waves production during inflation [ p^|j20| ]. The adiabatic theorem |l^,^ indeed implies that |/3fesP 
falls off exponentially with k beyond some cut-off kc (see also Ref. for a recent study). In effect, a numerical 
integration of the field equation shows that if the transition between inflation and radiation/matter domination is 
sufficiently smooth, an exponential cut-off appears, as shown in Fig. (^. Since the integral in Eq. ( pT| ) is proportional 
to kc (provided kc\'qi\ > 1), we will use numerical integration of the field equation for quantitative estimates, and the 
analytical solution to understand the behavior in various regimes of fc, /zi, and fj,2- 

In Fig. (|l|), we compare the analytical (dashed line) and numerical solutions to \(3ks\^ in the case where ^1 = 
and /i2 = 1, and for a transition into matter domination. The exact value of the cut-off wavenumber kc depends 
on the duration A77 of the transition between infiation and matter domination ||l^,0, as indeed, the "degree of 
non-adiabaticity" of the transition is inversely proportional to A77. Figure (|]) shows that the analytical solution is 
an excellent approximation to the numerical solution for k ^ kc, even for fcc|?/i| > 1. Furthermore, as expected, 
kc oc 1/Ar/, and therefore the number density in Eq. ( p^ also scales approximately as l/Ary. Indeed, as shown in 
Fig. (|l|), the analytical solution for the sudden transition corresponds to the numerical solution with the cut-off kc 
cast to infinity. 

In Fig. (H), we show a (base 10) logarithmic contour plot of Yrj^2i9 *s/ 200) (T^/Hi)^ {an /ai)'^ in the plane ^i,/i2, 
assuming a transition with Ar^/jTyil = 1 into matter domination. A transition into radiation domination gives similar 
results. As fJ-i = fJ-2 ^ +00, the production is exponentially suppressed [see also Fig. (H)], in agreement with the 
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FIG. 2. Logarithmic (base 10) contour plot of Yzji = y3/2(<7*s/200)(T'R/JJ"i)^(aR,/ai)^, for Ari/\i]i\ = 1, and for a transition 
into matter domination, in the plane pi — j-i2- 



adiabatic theorem. This can be seen in Eq. (|T6|), at least in the limit where fc|?7i| <C 1, for which |/3fesp oc exp(— 27r/ii). 
In the limit /ui — > 0, with ^2 fixed, one has ¥3/2 oc /x^, and similarly for ^2 — > 0, with /ii fixed. Notably, ^3/2 —> as 
/ii — > and ^2 0, since a massless gravitino is conformally invariant. Finally, for a fixed /ii, the number density 
is not exponentially suppressed as ^2 +00; indeed, Eq. (|l^) gives |/3fcsp ^1/2 when ^ and /i2 ^ +00, for 
A:|77i| < 1. Since kc is roughly proportional to max(/ii, ^2), the number density increases linearly with fj,2 in this limit. 
Note that this does not contradict the adiabatic theorem: as /.ti — *■ and fi2 +00, for a fixed A77, the transition 
becomes increasingly non-adiabatic, and particle production is not suppressed. 

Finally, in Fig. (^), we show cuts of the previous contour plot, for /ii — ^2 [diagonal of Fig. (^)], for — [y— axis 
of Fig. (||)], and for fJ,2 = [x— axis of Fig. (||)], in each case for a transition into matter domination with Ar//|?7i| = 1. 
For the sake of completeness, we also include the result of a numerical integration for the case where the mass of the 
gravitino is constant and has the same value in both inflationary and matter dominated epochs, in dash-dotted line. 
Just like the case /ii = /i2, it shows exponential suppression as m 3> -ffi. 

The number density of gravitationally produced gravitinos present at the time of reheating thus depends on several 
parameters, notably the effective mass terms during and after inflation, the duration of the transition from inflation 
to reheating, and the number of e-folds of reheating. The yield of gravitinos scales as the inverse of the transition 
timescale, which defines the degree of "non-adiabaticity" , and decreases as exp(— SiV^), where iV^ = ln(ar{/ai) is the 
number of e-folds of reheating, during which the gravitinos are diluted. These dependences make a direct comparison 
with the number density of gravitinos produced in particle interactions in reheating slightly delicate. 

Let us first isolate the dependence on the mass terms and transition timescale in the fiducial quantity 13/2, which 
is defined through: Y3/2 = y3/2(5*s/200)"i(i/i/rR)3 exp(-3iVR); the quantity plotted in Figs (|) and (|) is ^3/2 
(for A77 = H^^). The number of e-folds of reheating, and therefore I3/2, depend on the detailed mechanism of 
reheating, which is unfortunately not well known at present. In the most standard model of reheating |2^ , in which 
the infiaton slowly decays through its coherent oscillations, and the Universe is matter dominated, one obtains: 
3iVR ~ 58.5 In [(ffi/10i3QgV)2(rj^/109GeV)-4(5*/200)-i]. The dilution is therefore quite strong, and ^3/2 ~ 
3x 10~"'^'*y"3/2(^^i/10"'^^ GeV)(TR/10^ GeV). More generally, if reheating takes place in an era dominated by an equation 
of state of the form p = wp, the above number of e-folds is reduced by a factor 1/(1 -I- w); for w — 1/3, for instance, 
which corresponds to a relativistic fluid, one finds: ^3/2 ~ 9 x 10-^^3/2(^1/10" GeV)3/2(g(^^/200)-i/^_^and the 
number density produced becomes independent of the reheating temperature. It was pointed out in Ref. p9| that in 
a general case, oscillations of an inflaton in a potential ^ Ai/)^" would yield an equation state with w ~ (n — l)/(n + 1) 
after averaging out over an oscillation period. Thus, in particular, for chaotic inflation with a potential A0^, the 
Universe is indeed dominated by a relativistic fluid during reheating {w = 1/3). 

The ratio of the number density of gravitinos produced in particle interactions during reheating to the entropy 
density is, up to logarithmic corrections [0: 13^2 ~ 3.7 x 10""(Tr/10^ GeV)(g*s/200)"3/2^ Therefore, the ratio 



6 




0.01 0.10 1.00 10.00 

m/H 



FIG. 3. Dependence of Y3/2 = ¥3/2(9 *s/200){Tn/Hif [an/ aif on the effective mass of the gravitino in different cases: solid 
line m/H = fj,i = 112, dotted line m/H = /i2 and f^i = 0, dashed line m/H = fj,i and /i2 = 0, and in dash-dotted line, as a 
function of a constant mass term, equal in both inflation and matter domination eras. In each case, the transition is operated 
on a timescale Ar;/|r;i| = 1 into matter domination. 



^3/2/^3/2 of these two yields, assuming that reheating takes place in a matter dominated era is: 



^-o.iy,J^)m"\ (18, 



3/2/ -^3/2 

and according to Fig. (^, ~ 10^^^^ 2(^'7/l'7il)^^j if ^1 = and fi2 = ^ or the reverse. This constitutes our 
main result. If throughout reheating, the Universe is dominated by a relativistic equation of state, it becomes: 

i^"'"^'^^='/<ioT^j (lo^j (200) ■ 

Whereas the ratio of the two production yields is independent of the reheating temperature when the Universe 
is matter dominated during reheating, it becomes inversely proportional to Tr when w = 1/3 (relativistic fluid). 
Therefore, if gravitational production is efficient, i.e., if m ~ during or after inflation, a low reheating temperature 
does not exclude a strong gravitino problem. 

Let us now discuss the magnitude oi As seen in Fig. (|), one probably has 1^3/2 < lO^^ if At; = , 

where the upper limit corresponds to /xi = and /X2 = 1, or = 1 and /i2 = 0. Therefore, for reheating in a 
matter dominated era, one finds that the production of gravitinos out of the vacuum is less efficient than that during 
reheating, provided A77 ^ 10~^-ffj~^. In the other limit, where the Universe is dominated by a relativistic equation of 
state during reheating, one finds that gravitational production can be much more efficient that reheating production 
of gravitinos, by a factor - 10^^? 2(^'7/lml)"^(TRyiO^ GeV)~i. 

Finally, an order of magnitude estimate for A77 is 4)/4> (0 is the inflaton field) taken at the point at which the 
slow-roll approximation breaks down, i.e. where 0^/2 ~ y{4>) (a dot denotes differentiation with respect to cosmic 
time). This gives At; « 2(0/A/pi)|t7i|, with Afpi = (87r)^/^Tn,pi the Planck mass. Therefore, ATy ^ |t7i| for scenarios of 
the chaotic type, and At; <C_[t;i| for scenarios of the new inflation type with small fleld values; quite possibly, in this 
latter case. At; < 10-''|t7i| |20|. 

In new inflation, therefore, the gravitational production cannot be neglected if ;xi and ^2 ^ ^ (or the reverse), 
(f) < 10~'*Mpi at the end of slow-roll, and Hi ~ 10^^ GeV. Moreover, if reheating proceeds faster than in the "standard" 
model (matter domination), such as in chaotic inflation, gravitational production of gravitinos can become more 
efficient than reheating production. If gravitational production dominates, cosmological bounds on the gravitino 
abundance at the time of Big Bang nucleosynthesis should be turned into upper limits on the effective mass terms of 
the gravitino during and after inflation, as gravitational production is suppressed as [m/HiY if m < Hi. 
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IV. DISCUSSION 



We discussed the conformal behavior of the gravitino in a spatiaUy flat FRW background spacetime. We obtained 
the Unearized field equation for the helicity ±3/2 components, which reduces to a Dirac like equation in curved 
spacetime. A massive gravitino is not conformally invariant, and cosmological particle production ensues, through the 
amplification of the vacuum fluctuations by the non-static background metric. We assumed that the gravitino effective 
mass is proportional to the Hubble constant, and used the technique of Bogoliubov transforms to calculate the ratio 
^3/2 of the number density of gravitinos to the entropy density at the time of reheating. This quantity depends on the 
effective mass of the gravitino during and after inflation, on the Hubble constant at the exit of inflation {Hi), on the 
duration of the transition between inflation and radiation/matter domination (A77), and on the number of e-folds of 
reheating. Notably, scales as the inverse of the transition timescale, which defines the degree of "non-adiabaticity" 
of the transition during which gravitinos are produced. The comparison of the gravitational production of gravitinos 
to production in reheating depends on the details of the mechanism of reheating, during which the gravitationally 
produced gravitinos are strongly diluted. 

If we assume that the gravitino mass is of order of the Hubble constant during or after inflation, that Hi ~ 10^^ GeV, 
and that the Universe is matter dominated throughout reheating, gravitational production is generically less efficient 
than production in reheating interactions, provided the transition is not too abrupt, i.e. A77 ^ 1Q~^H^^. However, 
in scenarios of new inflation, one can find A77 < lO-^iffS in which case gravitational production would turn out to 
produce as many gravitinos, or more, than reheating interactions. Similarly, if reheating proceeds faster, for instance if 
the Universe is dominated by a relativistic fluid during reheating, as happens in e.g., chaotic inflation with a potential 
A(/>^, the number density of gravitinos produced out of the vacuum exceeds, possibly by a large factor, the density 
of gravitinos produced in particle interactions during reheating. It must be stressed that in the above, we assumed 
the effective mass m of the gravitino to be of order the Hubble constant, either during or after inflation. If not, 
gravitational production is suppressed as {m/Hif'. 

To conclude, let us briefly address the particular case of pre-Big-Bang string cosmology in which particle 
production out of the vacuum has been studied extensively [^,p4j, but not for spin— 3/2. In this scenario, one expects 
that, to leading order, the gravitino mass term vanishes during inflation, if the only dynamical field is the axion-dilaton 
field, as indeed, the tree level superpotential of string-inspired supergravity does not receive contributions from the 
dilaton. Similarly, m = also in the post-inflationary phase (if only the axion-dilaton field is considered), at least 
until a non perturbative superpotential for the dilaton sets in, or until supersymmetry breaking takes place. If we 
assume that the gravitino is also massless during the so-called stringy phase, it then couples to the axion-dilaton 
field only through the Kahler connection. The field equation Eq. (||) then decouples into two first order differential 
equations for the mode functions hss and gss^ whose solutions are written as in flat-space up to a time-dependent 
phase which depends on the Kahler connection. If the initial state corresponds to the conformal vacuum as — > — cx), 
then h™j^ = g"^_ = 0. Since we assume that the gravitino remains massless after inflation, conformal triviality also 
holds as ?7 +00, and h°^^ = g'^}, = 0. From Eq. (p^), it is then obvious that Pks = 0, i.e. no particle production 
takes place. 

At the next level of approximation, one should consider moduli fields, take into account higher order corrections to 
the effective action in the stringy phase, and/or introduce a non-perturbative superpotential to stabilize the dilaton 
in the FRW phase. This would lead, quite presumably, to the appearance of an effective mass for the gravitino. 
Unfortunately, these effects are difficult to implement, because the underlying dynamics or the physics remain poorly 
known. To give an example of what could be obtained, let us assume that the gravitino is massless during the pre-Big 
Bang and stringy phases, and that it acquires a mass after the exit in the FRW era. Then the methods and results 
of the previous section can easily be transposed to this scenario, since the gravitino, being massless in the pre-FRW 
eras, is insensitive to the background dynamics. It is easy to verify that if the exit in the FRW phase takes place 
at a scale Hi ~ 10^^ GeV, as has been advocated recently [Q, gravitational and reheating production of gravitinos 
become of the same order, even if the Universe is matter dominated throughout reheating, provided m ^ Hi. 

However, reheating in pre-Big Bang cosmology is not expected to proceed through coherent oscillations of the 
"inflaton" , and the above estimate could turn out to be naive. A detailed study of the mechanism of reheating in 
pre-Big Bang cosmology thus appears mandatory. Depending on how fast reheating proceeds, and what temperature 
is achieved, this could lead to a strong gravitino problem (which one would naively expect if Hi ~ 10^^ GeV), which 
would thus require: TO3/2 ^ 10'* GeV for an unstable gravitino, or 7713/2 ^ 2 keV for a stable gravitino. A more detailed 
study of this problem is left for further work. 

Note added 

Upon completion of this paper, we became aware of a related work by A. L. Maroto and A. Mazumdar ("Production 
of spin 3/2 particles from vacuum fluctuations", hep-ph/9904206D . These authors obtained the field equation for 
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helicity 3/2 gravitinos, assuming 7'*^'^ = (which projects out helicity 1/2 modes), and calculated the amplification 
of vacuum fluctuations, using the technique of Bogoliubov transforms. They applied their technique to the production 
of gravitinos in preheating. In this respect, their work and ours are complementary: gravitational production during 
inflation generically produces particles with wavenumber k Hi, while in preheating, the production takes place for 
modes with Hi. 

After the present paper was submitted, two other related studies appeared: R. Kallosh, L. Kofman, A. Linde and 



A. Van Proeyen ("Gravitino production after inflation", hep-ph/9907124) studied the problem of gravitino production 



during inflation and during preheating, for both helicity 1/2 and helicity 3/2 modes. Their important work shows 
that the helicity 1/2 modes are not conformally invariant even if they are massless, and that their production in 
preheating can be very large. The paper by G. F. Giudice, A. Riotto and I. Tkachev ("Non-thermal production of 



dangerous relics in the early Universe", hep-ph/9907510) reaches similar conclusions 
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APPENDIX: NOTATIONS 

We write 7^ a general relativistic Dirac matrix, and 7^ = 6^7^, or, if confusion could arise, 7a, a (constant) 
flat-space Dirac matrix. We deflne: aab — 5(70, 7h]- The Dirac matrices are written in the Weyl representation: 

with: (Jq — (1,<t), and CTq = (1,— cr), and the a are flat-space Pauli matrices. We also define: 75 = 47071 7273- 

Our choice of vierbein for the FRW background is: e° = 0(77), — a{r})~^ . The spin connection, without 
torsion, is then: 

i^o"Vafc = 0, ic^faab = 1^7.7°, (A2) 

where we defined H = a' /a. We define the helicity operator e^cr for a Weyl spinor of momentum fc, where e}" is the 
unitary vector along k. We then define X-i-(^) ^-nd X-(^) the eigenspinors of e^cr with respective helicity -1-1/2 
and —1/2: e^crx±(fc) = ±x±(fc). We decompose a four component spinor ^{r],k) in eigenstates of helicity [ pT| , 

-0(?7, k) = tp+ir], k) + tp-it], k), with: 

where hsirf), and gs(?7) are (scalar) functions of conformal time. The eigenspinors x+ and X- verify, in particular: 
x\xs' = ^ss' ■ In Section II, we perform a similar decomposition for the spinor-vector in terms of the mode functions 
hmsiil)i where m = L, +, — denotes the helicity of the polarization vector, and s = ± denotes the spinor helicity. 

Finally, we define the charge conjugation operator: C = 17^7", and the conjugate -0*^ of a spinor ip: ip'~^ = iC7°\l/*. 
One can show that: icf2X*s = ~sXs, where s = ±. It is then easy to show that the conjugate of a spinor ips{r],—k), 
with helicity s, and momentum — fe, is: 

^?{v,-k)^^^Jl^l). (A4) 
This identity is useful in deriving the normalization identities of the gravitino operator in Section II. 
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